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I . INTRODUCTION 

Computer  programs  for  R-field,  H-field  and  combined  field  solutions 
for  a perfectly  conducting  body  of  revolution  excited  by  an  obliquely 
incident  plane  wave  electromagnetic  field  are  described  and  listed  in  this 
report.  The  H-field  solution  is  obtained  by  applying  the  method  of  moments 
to  the  H-field  integral  equation.  The  E-field  solution  is  obtained  by 
applying  the  method  of  moments  to  the  E-field  integral  equation.  The  matrix 
equation  for  the  combined  field  solution  is  a linear  combination  of  the 
matrix  equations  for  the  H-field  and  E-field  solutions.  The  general  theory 
and  the  method  of  computation  are  given  in  reference  (13-  Equations  numbers 
drawn  from  reference  [1]  are  preceded  by  1-.  For  instance,  (1-J)  denotes 
equation  three  of  reference  [1]. 

The  computer  program  subroutine  YMAT  calculates  the  square  moment 
matrix  foi  the  H-field  solution  and  is  described  and  listed  in  Section  II. 

The  subroutine  ZMAT  calculates  the  moment  matrix  for  the  E-field  solution 
and  is  described  and  listed  in  Section  III.  The  subroutine  YZ  calculates 
the  moment  matrices  for  both  the  H-field  and  E-field  solutions  simultaneously 
and  is  listed  in  Section  IV.  The  subroutine  PLANE  calculates  the  plane  wave 
measurement  vectors  and  is  described  and  listed  in  Section  V.  As  shown  in 
[1]»  the  plane  wave  excitation  column  vectors  for  the  H-field  and  E-field 
solutions  and  hence  also  for  the  combined  field  solution  can  be  expressed  in 
terms  of  certain  plane  wave  measurement  vectors.  A main  program  which  uses 
YMAT,  ZMAT,  YZ,  and  PLANE  to  obtain  the  H-field  solution,  the  E-field  solu- 
tion, and  the  combined  field  solution  is  described  and  listed  in  Section  VII. 
The  final  section  shows  some  examples  of  computations  made  with  the  program 
and  gives  a discussion  of  them. 


[1]  J.  R.  Mautz  and  R.  F.  Harrington,  "H-Field,  E-Field  and  Combined 
Field  Solutions  for  Bodies  of  Revolution,"  Interim  Technical 
Report  RADC-TR-109,  Rome  Air  Development  Center,  Griff iss  Air 
Force  Base,  New  York,  March  1977 


II.  THE  SUBROUTINE  YMAT 


A.  Description; 

The  subroutine  YMAT(NN,NP,NPHI,RH,ZH,X,A,Y)  stores  the  matrix 


[Y] 


,tt 

Yt* 

n 

n 

4 1 

yM 

n 

n 

(1) 


appearing  in  the  H-field  matrix  equation  (1-17)  by  columns  in  Y.  The 
elements  of  [Y]  are  given  by  (1-31).  The  input  variables  are  defined 
in  terms  of  variables  appearing  in  reference  [1]  by 


NN  - n 
NP  = P 
NPHI  - N 

<t> 

RH(i)  - kp~ 
ZH(i)  = kz~ 
X(k)  *= 

A(k)  **  A^ 


In  summary,  n denotes  e^n^  dependence,  (p^,  z^) , i = 1,2,...P  are  coordi- 
nates on  the  generating  curve,  the  k which  multiplies  and  z^  is  the 
propagation  constant,  and  x^  and  A^  are  the  abscissas  and  weights  for  the 
N point  Gaussian  quadrature  integration  in  $. 


Minimum  allocations  are  given  by 


COMPLEX  Y (4*N*N) 

UIMENSION  RH(NP),  ZH(NP),  X(NPHI) , A(NPHI) , 
D(NG),  PD(NG),  CR(NPHI),  Cl (NPHI), 
C2(NPHI),  C3(NPHI) 

COMMON  RS(NG) , ZS(NG),  SV(NG),  CV(NG) , T(4*N) 
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where 


NG  = NP  - 1 
N = (NG-2)/2 

The  variables  In  common  make  the  results  of  some  intermediate  calculations 
done  in  YMAT  available  to  the  subroutine  PLANE  described  in  Section  V,, 

no  loop  57  puts  kd^  of  (1-28)  in  n,  kp^  of  (1-27)  in  RS,  kz_^  of  (1-27) 
in  ZS,  sin  v^  in  SV,  cos  v^  in  CV,  and  - in  PD.  DO  loop  68  puts  (1-29)  in  T. 


^ipi 


2 A, 


With  regard  to  (1-11)  and  (1-24)  - (1-26),  DO  loop  25  puts  4 sin  (-—) , 

nA^Bin2 (~~)  cos(n<t>k),  j A^cos  ^ cos(n<j>k),  and  — A^  sin  ^ sinfn^)  in  CR, 
Cl,  C2,  and  C3  respectively  where,  as  prescribed  by  (1-37), 


*k  ' 2 (\  + » 

DO  loop  62  initializes  Y. 

The  computation  of  (1-31)  is  sequential  in  s rather  than  i) . A 
computation  sequential  in  s is  preferable  because  it  does  not  require 
storing  (Yl)  , (Y2)  , (Y3)  , or  (Y4)  versus  s.  Unfortunately,  (1-31) 
as  it  stands  tacitly  implies  computation  sequential  in  ij . An  alternative  form 
of  (1-31)  somewhat  more  suggestive  of  computation  sequential  in  s is  given  by 


1 

✓—N 

rt 

rt 

w 

1 

n (Yl), 

n pq 

ij 

(Y*1) 

o pq 

2p+2  2q+2 

(Y2)i1 

<yA 

" . J T2p-2+i  J T2q-2+i 

i“2p-l  r j=2q-l  M ' 

(Y3),, 

r-  pq 

i.1 

(VM) 

(Y4) 

n pq 

L 

(2) 


If  (1-38)  is  true,  then 

((Y“) 
n pq  i.1 

((vf)  )., 

n pq  i.1 


- ((Y**)  ) 

u n qp  ji 

- ((Y*C)  ),. 

n op  1i 


(Or**)  ).,  = - ((rt<p)  ) 

n pq  ij  n qp  ji 


(3) 
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where  the  additional  subscripts  ij  denote  the  contribution  due  to  the 
i,1tk  term  on  the  right-hand  side  of  (2).  From  (3) 


(Y“> 

«r> 

(Y**) 

n 


pq 


pq 


pq 


pq 


((Y“)nn),<, 

n pq  i£j 


((Y 

n pq  i<j 


U n 


U n ;pq;isl 


u n ;qp;i<j 


((Yf)aD)i<, 
n qp  i£j 


((Y^)  ),<, 

n qp  i<j 


((Ytt)  ) 

n ;qP;i<J 


(4) 


where  the  subscript  notation  i<j  denotes  the  sum  of  the  contributions 
due  to  all  terms  for  which  i <•  j plus  half  of  the  i = j term  on  the 
right-hand  side  of  (2).  According  to  (4),  it  is  sufficient  to  calculate 
only  the  i<j  terms  on  the  right-hand  side  of  (2) . Inspection  of  the 
limits  of  summation  in  (2)  shows  that  the  first  term  on  the  right-hand 
side  of  (4)  is  zero  if  p > q+1  and  that  the  second  term  is  zero  if 

p < q-i. 

The  ij*"  term  of  the  first  of  equations  (2)  contributes 

tt  k=1’2 

T2(p+k-l)+i  12(q+£-l)+1(Y1)ij  t0  (Yn  } (p+k)  (q+£)  * (5 

where  p « [(i+l)/2]-2 

q = l(j+l)/2]-2 

where  [1]  denotes  the  largest  integer  which  does  not  exceed  i.  The  k=l 
term  should  be  omitted  from  (5)  for  the  first  two  values  of  i,  and  the 
k=2  term  should  be  omitted  for  the  last  two  values  of  i because  the 
triangle  functions  do  not  overlap  at  the  ends  of  the  generating  curve. 
Similarly,  the  £-1  term  should  be  omitted  for  the  first  two  values  of  j, 
and  the  1^2  term  should  be  omitted  for  the  last  two  values  of  j.  The 
last  three  of  equations  (2)  imply  relations  similar  to  (5) . 
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The  indices  I and  J of  nested  DO  loops  60  and  59  are  respectively 

i and  j of  (5).  DO  loop  61  accumulates  G^,  , and  G^  of  (l-24)-(l-26) 

for  I )*  J and  . 5G^,  . 5G2*  and  .50^  for  I * J in  Gl,  G2,  and  G3.  The 

(Yl)^,  (Y2)ij,  (YS)^  and  (Y4)^  which  appear  in  (5)  and  similar  equa- 
tions and  which  are  given  by  (1—32)  — ( 1—  35)  with  the  r- terms  missing 

kdipi 

and  the  i*j  result  halved  are  put  in  Yl,  Y2,  Y3,  and  Y4  just  after  DO 
loop  61.  The  indices  K and  L of  nested  DO  loops  32  and  31  are  respec- 
tively k and  £ of  (5).  The  Y(J1),  Y(J2),  Y(J3),  and  Y(J4)  in  DO  loop 

32  are  respectively  «“)  (p4k)  (q+t)  . Y*b  (p+k)  (,+l)  . «**>  (ptk)  (,+ll , . 


and 


(p+k)  (q+O ' 


DO  loop  11  carries  out  (4)  for 


( (J,J)  J = 1 

<p,q)  “ \ 

/ (J,J),  (J-1,J),  and  (J,J-1)  J > 1 


DO  loop  11  also  adds  the  contributions  to  (2)  due  to  the  t-tt terms 

k diPi 

in  (1-32)  and  (1-35).  In  DO  loop  11,  Y(KD1),  Y(KD2),  Y (KD3) , and 

Y(KD4)  are  the  diagonal  (J,J)  elements  of  Ytt:,  Y1^,  Y**  and  Y^^  respec- 

n n n n 

tively.  Similarly,  Y(KU1),  Y(KU2),  Y(KU3),  and  Y(KU4)  are  the  (J-1,J) 
elements  while  Y(KL1),  Y(KL2),  Y(KL3) , and  Y (KL4)  are  the  (J,  J-l)  ele- 
ments. The  first  R1  in  DO  loop  11  is  the  contribution  to  the  diagonal 

elements  of  (2)  due  to  the  ^7-^--  terms.  The  second  R1  in  DO  loop  11  is 

the  contribution  to  the  off  diagonal  elements  due  to  the  r-7-2 terms. 

k dipi 

Nested  DO  loops  13  and  14  use  (1-39)  for  (i,j)  = (I,J)  to  put 

0?“)^,  (YJt)i y and  >Vn<|))ij  for  1 i1'2  in  Y<KL1>»  Y(KL2), 

Y(KL3),  and  Y(KL4)  respectively. 
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B 


LISTING  OF  THE  SUBROUTINE  YMAT 


SUBROUTINE  Y H A I ( NN ,N P , NPHI , RH , ZH .X , A, Y ) 

COMPLEX  U*Y(1600)»GlfG2»G3fYl»Y2fY3*Y4 

DIMENSION  PH (A3) ,ZH(A3) ,X(20) , A ( 20 ) , D( A2 ) ♦ PO ( A2 ) ,CP (20) tC 1(20) 
OIMENSICN  C 2 ( 20) tC3( 20 ) 

COMMON  RS(A2),ZS(A2) , SV ( A2 ) t CV ( A2 ) t T ( 80 > 

Pl=3.1A15S3 


DO  57  1 = 2, NP 

12=1-1 

DR=RH( I )-RH  ( 12) 

DZ=ZH( ( )-ZH( 12) 

0(12  )"SQRT(DR*DR-*-DZ*DZ) 
RSU2)=.5*(RH<  I)*RH(  12)) 
ZSU2)  = .5*(ZH(  I)*ZH(  12)  ) 
SV( I 2) =0P/0 ( 12) 

CV( I2)=CZ/C( 12) 
PC(I2)=PI/(C(I2)*RS( 12)) 

57  CONTINUE 
NG=N  F-l 
N2=NC— 2 
N=N2/2 
Jl  = l 
J5=l 

DO  68  J=  1, N 
D1=C(J1) 

D2=D  < J 1 ♦ 1 1 
C3=C  I J 1 *2 ) 

DA=D ( J1 +3 ) 

DEL l=Dl*D2 

0EL2  =03-fCA 

T( ) = • 5*01*0  1 /DELI 

T(J5*1 )=(D1+,5*D2)*02/0EL1 

T( J5*2)*(0A*.5*C3)*03/0EL2 

T(J5+3I=.5*CA*CA/DEL2 

ji=ju? 

J5=J5+<t 
68  CONTINUE 
P 12=  ,5*PI 
FN  = NN 


DO  25  K= 1 1 N PH  I 
PH=P  12* (X ( K I + 1 • ) 

PHN=PH*FN 
SN=S  IN  ( .5*PH  ) 
CR(K)=A.*SN*SN 
R l=P 12  *A ( K > 

Cl(K)=.5*Rl*CP(K) *COS( PHN ) 
C2(K)=R1*C0S(FH)*CCS(PHN) 
C3(K  )=Rl*SlN(Ph)*SIN(PHN) 
25  CONTINUE 
N2N=N2*N 
NAN=N2N*2 
DO  6?  J= 1 , N AN 
Yl J)=0. 

62  CONTINUE 
U=  ( 0 . , 1 . ) 

CO  55  J= 1 » NC 

LI  = 1 

L2=2 

IFU.LE.21  Ll=2 
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IFU.GT.N2)  12  = 1 
J8  = J + l 
J9=J8/2 
JT=2*J9+J~6 
J5  = < J9-3  )*N2-2 
Sl=l. 

DO  60  1=1, J 
18=1+1 
19=18/2 
IT=2*l9+I-6 
J6=I 9+  J5 
RP=RS l J )-RS( I ) 

ZP=ZS( J)-ZS(I) 

R2=RP*RP+ZP*ZP 

IF ( I .NE . J ) GO  TO  41 

Sl  = . 5 

R2=. C625*D( J ) *C ( J) 

41  R3=RS(IJ*RS(J> 

G 1=0. 

G2=0. 

G3=0  . 

00  61  K=  l , NPH I 
R4=R2+R3*CR(K) 

R5=SCRT(P4) 

Yl=Si/(R4*R5i*<  l.+U*R5)*ICQSIR5)-U*SIN<R5» > 

G1=CIIK)*YI+G1 
G2=C2( K J + Y1+G2 
G3=C3( K ) *Y 1 +G3 
61  CONTINUE 
G3=U*G3 

Y1=(PP*CVU)-ZP*SV(J))*G2-RS  !I»*CV(J)*Gl 

Y2  = I RS(  J)*SVU)*CVU)-RSm*SV(J)*CVll  )-ZP*SV(I)*SV(  JJ)*G3 
Y3=ZF*G3 

Y4= ( RP*CV ( I J-ZP+SV II )>*G2+RS  f J)*CV( I >*Gl 
Kl=l  f 

K 2=2 

IFU.LF.2)  Kl  = 2 
IFII.GT.N2)  K2  = l 
DO  31  L=L1,L2 
LT=JT+L+L 
J7= J6+L  *N2 
00  32  K=Kl , K2 
KT=I T+K+K 
TT=T  (LT  l*T ( KT  ) 

J1=J7+K 
J2= J 1+N 
J3= J 1+N2N 
J4=  J3+N 

Y ( J 1 >=  T T +Y  1 + Y I J l ) 

Y( J2 )=IT*Y2+Y( J2  ) 

Y(  J 3 )=T  T *Y  3 ♦ Y ( J3) 

Y(J4  )=TT  +Y4+Y  ( J4  ) 

32  CONTINUE 
31  CONTINUE 
60  CONTINUE 
59  CONTINUE 
N 2 P=  N 2 + 1 
KDl  = l 
J 1=  1 
J5  = l 
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DO  11  J=1,N 
KD2=KD1+N 
KC3=KOl+N2N 
KC4=KD3+N 
R2=T  I J 1 1 » 

R3=T ( J 1 + 2) 

R4=T  (Jl  + 3) 

J6= J5+ 1 

R 1=T  ( J 1 ) *T(  J 1 ) *PD ( J5)+R2*R2*PD( J6)+R3*R3*PD( J6+1  ) +R4*R4*PC { J6+2 ) 
G1=Y (KDi )-Y  ( KD4  ) 

Y ( KD l I =R l+G 1 
Y(KQ2)=0. 

Y ( KD3 ) =0  • 

Y(  KD4)  = R 1— G 1 

IFU.EQ.U  GO  TO  22 

KUl=KDt-l 

KU2=KD2—  1 

KU3=KD3-l 

KU4=KD4~1 

KLl=KDl-N2 

KL2=KD2-N2 

KL3=KD3~N2 

KL4=  KD4-N2 

Rl=T( J1I*T( Ji-2)*PD( J5)+R2*T (J1-1)*P0< J6) 

GI=Y(KU1I-Y(KL4  ) 

G2=Y  ( KU4 )-Y ( KL  1 ) 

Y(KU1)  =R  l+G  1 
Y(KU2)-Y(KU2)-Y(KL2) 

Y(KU3l=Y(KU3)-Y(KL3) 

Y( KU4)=R 1+G2 
Y(KL1)=R1-G2 
Y ( KL21 =— Y( KL2I 
Y ( KL  3 )=-Y  t KL3 ) 

Y ( KL4I =R1-G1 
22  KDl=KDl+N2P 
Jl=J 1+4 
J5=J  5+2 
11  CONTINUE 

IFCN.LT.3)  RETURN 
J2=N2 

00  13  1=3, N 
J2= J2+N2 
J 1=1-2 
KL1=  I 

00  14  J=l,Jl 
KUi= J2+J 
KU2=KU1 +N 
KU3=KU 1+N2N 
KU4= KU3+N 
K12=KL1+N 
KL3=KL 1+N2N 
K14=KL3+N 
Y (KL l ) =— V ( KL4  ) 

Y(KL2)=-Y( KU2  ) 

Y ( KL3 ) = -Y ( KU3  ) 

Y(KL4) =-Y (Kill 
KL 1=  KL 1 +N2 
14  CONTINUE 
13  CONTINUE 
RETURN 
END 


a 


III.  THE  SUBROUTINE  ZMAT 


A.  Description 


The  subroutine  ZMAT  (NN,  NP,  NPHI,  RH,  ZH,  X,  A,  Z)  stores  the 


matrix 


[Z] 


NP, 

NPHI,  : 

,tt 

'n 

n 

M 

Z** 

'n 

n 

(6) 


appearing  in  the  E-field  matrix  equation  (1-57)  by  columns  in  Z.  The 
elements  of  [Z]  are  given  by  (1-69) -(1-72) . The  input  variables  are 
defined  in  terms  of  variables  appearing  in  reference  [1]  by 


NN  = n 
NP  =■  P 
NPHI  - N 
RF(i) 


* 

kp. 


ZH(i)  = kzi 
X(k)  = x^ 
A(k)  * A^ 


In  summary,  n denotes  e 


jn$ 


dependence,  (p^,  z^) , i=l,2,...P  are  coordi- 


nates on  the  generating  curve,  the  k which  multiplies  p^  and  is  the 

propagation  constant,  and  x^  and  A^  are  the  abscissas  and  weights  for  the 

N,  point  Gaussian  quadrature  integration  in 
9 

Minimum  allocations  are  given  by 
COMPLEX  Z(4*N*N) 

DIMENSION  RH(NP),  ZH(NP),  X(NPHI),  A(NPHI),  D(NG) , 

TP (4*N) , CR(NPHI),  C2 (NPHI) , C3(NPHI),  C4(NPHI) 

COMMON  RS(NG),  ZS(NG),  SV(NG),  CV(NG),  T(4*N) 
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NJ|=i 


where 


NG  = NP-1 
N = (NG-2)/2 

The  variables  in  common  make  the  results  of  some  intermediate  calculations 
done  in  ZMAT  available  to  the  subroutine  PLANE  described  in  Section  V. 

DO  loop  57  puts  kdj  of  (1-28)  in  D,  kp^  of  (1-27)  in  RS,  kz^  of  (1-27) 
in  ZS,  sin  v^^  in  SV,  and  cos  v.^  in  CV.  DO  loop  68  puts  (1-68)  in  TP  and 

2 ^k 

(1-29)  in  T.  With  regard  to  (l-62)-(l-65)  , DO  loop  25  puts  4 sin  (y) , 

cos  <(>k  cos  (n<J)k),  ^ A^  sin  sin  (n^k> , and  j Afc  cos  (n<j>k)  in  CR,  C2, 
C3,  and  f.4  respectively  where,  as  prescribed  by  (1-75), 

\ ■ f <\+1) 


DO  loop  62  initializes  Z. 

Equation  (1-69)  can  be  rewritten  as 


2p+2  2q+2 

>P,  ' 3 }X-l  'V**  WV,1#  ViSl"  + V“  ViC°'  V 


T t * n 1 

i2p-2+i  2q-2+j  4; 


(7) 


where  G.  and  G_  are  evaluated  at  (p,z,p',z’)  = (p . ,z. ,p. ,z  ) . From  (7)  and 
4 5 , . ..  1 1 J J 

similar  equations  for  (?,  ) , (Z  ?)  , and  (Zw)  , we  obtain 

n n pq  n pq  n pq 


(8) 


((Z**)  - («♦*)„„)., 

n pq  ij  n qp  ji 
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where  the  additional  subscripts  i,j  denote  the  contribution  due  to  the 
til 

ij  term  on  the  right-hand  side  of  (7)  and  similar  equations.  Equa- 
tion (8)  implies  that 


<c>p,  ■ «»£*>„>*„ + «zrvi<j 


(2+t) 
n pq 


n pq  ifj  n ' qp' i<j 


(9) 


(zf)M 

n pq 


n pq  i<j 


■ «zr>Pq>i<j + <(zr>qP>i<j 


where  the  subscript  notation  i<j  denotes  the  sum  of  the  contributions 
due  to  all  terms  for  which  i<j  plus  half  of  the  i=j  term  on  the  right- 
hand  side  of  (7)  and  similar  equations.  According  to  (9),  it  is  suf- 
ficient to  calculate  only  the  i<j  terms  on  the  right-hand  side  of  (7) 
and  similar  equations.  Inspection  of  the  .limits  of  summation  in  (7) 
and  similar  equations  shows  that  the  first  term  on  the  right-hand 
side  of  (9)  is  zero  if  p > q+1  and  that  the  second  term  is  zero  if 
p < q-1 . 

The  ij1"*1  term  of  (7)  and  similar  equations  contributes 


{T2(p+k-l)+iT2(q+2,-l)+j(Z1)  + T2(p+k-l)+iT2(q+il-l)+j(Z1A)}  t0  (Zn  ) (p+k)  (q+H)  ’ 

T2(p+k-l)+i{T2(q+£-l)+j(Z2)  + T2(q+i-l)+j(Z2A)}  t0  (Z^  ) (p+k)  (q+O  ’ 

t4  <*<>) 

T2(q+£-l)+j{T2(p+k-l)+irZ3)  + T2(p+k-l)+i(Z3A)}  CO  (Zn  } (p+k) (q+*) ’ 

flnd  44  1 k=1’2 

T2(p+k-l)+iT2(q+£-l)+.i(Z4)  t0  *"Zn  ^ (p+k)  (q+£)  . * |a=1>2 
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where 


Z1  = j(GjSin  v^sin  + G^cos  v^cos  v^) 
Z1A  = -jG. 


Z2  = - G.sin  v, 
6 j 

Z2A  = - ~ G. 

kpt  4 

Z3  * G,sin  v, 
o i 


Z3A  = G. 
kPj  4 


(ID 


24  " 3<S  - V 

p = 1(1+1) 12]  - 2 

q « ((.1+l)/2]  - 2 

where  [i]  denotes  the  largest  Integer  which  does  not  exceed  i.  The  k=l 
term  should  be  omitted  from  (10)  for  the  first  two  values  of  1,  and  the 
k=2  term  should  be  omitted  for  the  last  two  values  of  i because  the  tri- 
angle functions  do  not  overlap  at  the  ends  of  the  generating  curve. 
Similarly,  the  l = 1 term  should  be  omitted  for  the  first  two  values  of 
j,  and  the  1 = 2 term  should  be  omitted  for  the  last  two  values  of  j. 

The  indices  I and  J of  DO  loops  60  and  59  are  respectively  i and 
j of  (10).  DO  loop  61  accumulates  G . , Gc , and  G,  of  (1-62) -(1-64)  for 

H J O 

I J-  J and  .5G^,  ,5G^,  and  .5G^  for  I = J in  G4,  G5,  and  G6.  The  Z’s 
defined  by  (11)  for  1 4 J and  these  Z's  divided  by  two  for  T = J are 
put  in  Zl,  Gl,  Z2,  G2,  Z3,  G3,  and  Z4  just  after  DO  loop  61.  The 
indices  K and  L of  nested  DO  loops  32  and  31  are  respectively  k and 
2.  of  (10).  The  Z(J1),  Z(J2) , Z(J3),  and  Z(J4)  in  DO  loop  32  are 


In  (10), 
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respectively  (Ztt:),  .w  . , (Z^),  ,1W  , nN , (Z^),  .,  w and 

n (p+k)  (q+J.)  n (p+k)(q+£)  n (p+k)  (q+St)  * 

Un  '(p+k)(q+0* 


DO  loop  11  carries  out  (9)  for 


Cp*q) 


(J,J)  J = 1 


(J,J),(J-1,J),  and  (J,J-1)  J > 1 


In  DO  loop  11,  Z(KD1),  Z(KD2),  Z(KD3),  and  Z(KD4)  are  the  diagonal 
(J,J)  elements  of  Z^1",  Z^,  Z^,  and  respectively.  Similarly, 

Z(KU1),  Z(KU2) , Z(KU3),  and  Z(KU4)  are  the  (J-1,J)  elements  while 
Z(KLl),  Z(KL2) , Z(KL3),  and  Z(KL4)  are  the  (J,J-1)  elements. 

Nested  DO  loops  13  and  14  use  (1-76)  for  (i,j)  « (I,J)  to  put 

(Z^),.,  (***)..,  <Z*V,,  and  (Z**).A  for  i < i-2  in  Z(KL1),  Z(KL2), 
n i j n l j n i]  oil  — 

Z(KL3) , and  Z(KL4)  respectively. 


LISTING  OF  THE  SUBROUTINE  ZMAT 


SUBROUTINE  ZPA7(NN»NP»NPHI *RH*  ZH  »X ♦ A »Z ) 

COMPLEX  U,  Z{l6C0)tGltG2tG3tG4»G5fG6,Zl»Z2tZ3tZ<« 

DIMENSION  RH ( 43 ) fZHI43)fX(20)»A(20)tD(42)tTPI80)fCRI20)tC2(2C) 
DIMENSION  C 3 ( 2C ) *C4l  20) 

CCMMCN  RSI  42)  tZS(42) , S V ( 42 > , CV ( 42 ) t T 1 80 ) 


PI=3. 141593 
DO  5 7 I = 2tNP 


12=1-1 

DR=RH( I )-RH I I 2 ) 

CZ=ZH( I )-ZH(  12) 

0U2  )=SQRT(DR*DR*DZ*DZ) 
RS(  I2)=.5*tRH(I)«-RH(  12)) 
ZSU2)  = .5»(ZM(I)+ZH(  12)) 
SV( I 2) =DP/0( 12) 

CVU2)  = 0Z/DU2) 

57  CONTINUE 
NG=NP- l 
N2=NG-2 


N=N2 /2 


Jl=l 
J5  = l 

DO  68  J=1»N 
D1=D  U 1 1 
02=0  t J 1 ♦ 1 1 
D3=D ( J 1+2) 

C4=0  ( J 1 *3 l 
DELI =01*02 
DEL2=D3*04 
J6=J5*1 
J7= J6+  1 
J8- J 7*1 

TP( J5)=D1/0EL1 
TP ( J 6) =D2/DEL1 
TPIJ7)=~C3/CEL2 
TP(J8)=—  D4/C6L2 
T(J5)=.5*D1*TP(J5) 

T { J6  )»I01*.5*C2)*TPI J6) 

T(  J7)=-(04*.5*13)*TPIJ7) 
T{J8)=-.5*D4*TP(J8) 

Jl=J 1*2 
J5=J  5*4 
68  CONTINUE 
P I 2= • 5*P! 

FN=NN 

DO  25  K=l»NPHI 
PH*PI2*CXIK)*l.) 

PhN=PH*FN 
SN=S  INI  .5*PM 
CR(K  )=4.*SN*$N 
Rl  = P I2*AIK) 

C2iK)=Rl*C0SIPH)*C0SIPHNI 
C3IK  )=R1*SIN|PE)*SINIPHNI 
C4(K)=R1*CCS(PEN) 

2 5 CONTINUE 
N2N=N2*N 
N4N=N2N*2 
DO  62  J=  1 » N4N 
ZIJ)  =0. 
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6 2 CONTINUE 
U=<0.,1.) 

DO  5S  J= 1 » NG 
FJ=FN/RS (J ) 

L 1 =1 
L 2=2 

I F ( J ,L  E . 2 ) 11  = 2 
IFU.GT.N2)  1.2=1 

J8=j4i 

J9  = J8/2 
JT=2*J9* J— 6 
J5= ( J9-3  )*N2-2 
Sl=l  . 

CO  6C  1 = 1, J 

18=141 

19=1 8/2 

IT=2*I94l-6 

J6=I94J5 

FI=FN/RS(n 

RP=RSU)-RSU  ) 

ZP=ZS( J)-ZS(  I J 
R2=RP*RP+ZP*ZP 
I F I I .NE.J)  GC  TO  41 
Sl=.  5 

R2».Q625*C( 

41  R3=RS( l )*RS<J) 

G4=0  • 

G5=0. 

G6=0. 

DO  61  K*l,NPHf 
R4=R24R3*CR  IK) 

R 5=  SCR  T ( P4) 

Zl=Sl/R5#(COS(R5)— U*SIN(R5) ) 
G4=C4(K)*Z14G4 
G5=C2(K)4Z14G5 
G6=C3( K ) * 2 1 «G6 
61  CONTINUE 

Z1*U*(  SV(1)  *SVU)*G54CV(  I)*CV(  J)«G4) 

Gl=-U*G4 

Z2=—  SV( J l*G6 

G2=-Fl*G4 

Z3=SV(  U*G6 

G3=F  J*G4 

Z4=U*(G5-Fl*G3) 

Kl  = i 
K2=2 

!F(  I . L E • 2 ) K 1 = 2 
IFU.GT.K2)  K2  = i 
OH  31  L = L l , L 2 
LT=JT4L4L 
J7=J64L4N2 
DO  32  K=Ki,K? 

KT= 1 T4K*K 

TT  = T (LTHT(KT) 

J1  = J 7 ♦ K 
J2=J14N 
J3=J l 4N2N 
J4= J3+N 

ZU1  1=TT4Z14TP1LT)*TP1KT)*G14Z(  Jl) 

Z( J2  >=TT*Z24TF<LT)*T (KT ) *G2 4 Z t J 2 ) 


IS 


Z( J3  }=TT*Z3+TP( KT) *T 1 L T ) *G3+  Z ( J 3 1 
Z I J4  )=TT*Z4  + Z ( J4) 

32  CONTINUE 
31  CONTINUE 
60  CONTINUE 
5<3  CONTINUE 
N2P=N2*1 
KOl=l 

DO  11  J = 1,N 
Kf2=K01+N 
K03=KD 1+N2N 
K04=K03*N 

Z(K01)-Z(KD1)*2(K01) 

Z ( KD2) =Z ( KD2 )-Z { KD3 ) 

Z ( KD3 )=— Z ( K02  ) 
ZiKDAJ=ZIKD4)+Z(KD4) 

IF(J.EQ.l)  GO  TO  22 

KUi-KDl-1 

KU2=KD2-i 

KU3=K03-1 

KU4=KD4-1 

KL 1=  KD 1-N? 

KL2=KD2-N2 

KL3=KD3-N2 

KL4=KD4-N2 

Z(KUn=Z(KUlMZ(KLl) 

ZIKU2) =Z(KU2)-Z(KL3) 
Z(KU3)=ZIKU3)-Z(KL2) 

ZI KU4)=Z<KU4)*Z(KL4> 

Z(KLl>=ZiKUl) 

Z { KL2 ) =-Z ( KU  3 ) 

ZIKL31 =-ZIKL2) 

ZIKLA)=Z(KUA) 

22  KDl=K0l*N2P 
11  CONTINUE 

IFiN.LT.3l  PFTURN 
J2=N2 

on  13  1=3, N 
J2= J2+N2 
J 1=  I -2 
KL  1=  I 

00  l A J=i,Jl 
KU1=J2*J 
KU2=KU1*N 
KL3=KUl +N2N 
KUA=KU3*N 
KL2=KLl ♦N 
KL3=KL1+N2N 
KL4=  KL3+N 
Z(KLI)=Z(KU1 > 

ZiKL2)=-Z(KL3J 

ZiKL3)=-ZiKU2) 

Z ( KL  A I =Z I KUA  ) 

KL  i=KL  1+N2 
14  CONTINUE 
13  CONTINUE 
BFTUPN 
END 
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IV.  THE  SUBROUTINE  YZ 


A.  Description: 

The  subroutine  YZ(NN,  NP,  NPHI,  RH,  ZH,  X,  A,  Y,  Z)  stores  by 
columns  [Y]  of  (1)  appearing  in  the  H-field  matrix  equation  (1-17)  and 
{Z]  of  (6)  appearing  in  the  E-field  matrix  equation  (1-57)  in  Y and  Z 
respectively.  The  subroutine  YZ  is  the  subroutines  YMAT  and  ZMAT  com- 
bined into  one  subroutine.  If  both  [Y]  and  [Z]  are  required,  some 
compile  time  and  execution  time  can  be  saved  by  using  YZ  instead  of 
YMAT  and  ZMAT.  However,  YZ  requires  that  [Y]  and  [Z]  be  stored 
simultaneously  whereas  YMAT  and  ZMAT  do  not  necessarily  require  separate 
storage  locations  for  [Y]  and  [Z], 

The  input  variables  NN,  NP,  NPHI,  RH,  ZH,  X,  and  A have  the  same 
meaning  as  in  the  subroutines  YMAT  and  ZMAT. 

Minimum  allocations  are  given  by 

COMPLEX  Y(4*N*N) , Z(4*N*N) 

DIMENSION  RH(NP) , ZH(NP) , X(NPHI),  A(NPHI) , D(NG) , 

PD(NG) , TP (4*N) , CR(NPHI) , Cl (NPHI) , C2(NPHI), 

C3(NPHI) , C4(NPHI) 

COMMON  RS(NG),  ZS(NG),  SV(NG),  CV(NG) , T(4*N) 

where 

NG  =*  NP-1 
N = (NG-2)/2 

The  variables  in  common  make  the  results  of  some  intermediate  calculations 
done  in  YZ  available  to  the  subroutine  PLANE  described  in  Section  V. 

Because  all  variables  used  in  YZ  can  be  traced  back  to  YMAT  and 
ZMAT,  anyone  who  has  gone  through  YMAT  and  ZMAT  should  have  no  trouble 
with  YZ. 


B.  LISTING  CF  ThF  SUBROUTINE  YZ 

StPRCUT INE  Y2 ( NN» NP»  NPH I* RH» ZH,X,A,Yt Z ) 

CGMPLEXU,Y( 1600) ,Z( 1600) ,GlfG2,G3,G4,G5,G6,YltY2,Y3,Y4,Zl,Z2,Z3,Z4 
DIMENSION  RM43),ZH(43),X(20I,A(  20 1 , DC  42) , PD ( 42)  «TP(80)  ,CR(2  0I 
DIMENSION  Cl(20)tC2(20ltC3(20)fC4(20) 

COMMON  RS(42),ZS{42) , $ VC  421 , C V ( 4 2 ) ,T ( 80 ) 

PI=3 .141593 
DO  57  1=2, NP 
12=1-1 

DR  = RM  I )-RH(I2  ) 

DZ  = ZH( I l-ZH  C 121 

DC  12  )=SQRT(OR*CR+DZ*CZ  ) 

RS( I2)=»5*(RH(I )+RH( 12)1 
ZSU2)=.5*(ZH(I)*ZH(I2n 
SV ( 1 2 ) =CP/0  C 12 ) 

CVC I 2) =DZ/D (12) 

PD(I2)=PI/(C(12)*RS(  12)) 

57  CONTINUE 
NG=NP-  1 
N2=NC~2 
N=N2/2 
Jl=l 
J5  = L 

DO  68  J=1,N 
Dl=D ( J 1 ) 

D2=D( Jl+1) 

03=D( J 1+2 ) 

D4=D(J1*3) 

DEL1=0U02 
0EL2=03+C4 
J6=J5U 
J7=J6«- 1 
J8= J 7*1 

TP ( J5 ) =01 /DEL  1 

TP ( J6 ) =D2/DEL 1 

TPC J7)=-C3/CEL2 

TPCJ8)=-D4/CEL2 

TCJ5  )=.5*D1*TP(J5I 

T ( J6)  = (Di+.5*C2)*TPC  J6» 

T(J7)=-(D4+.5*D3I*TP(J7) 

TCJ8 )=-.5*04*TP( J8) 

Jl  = J 1+2 
J5= J 5+4 
68  CONTINUE 
P I 2= . 5*P I 
FN=NN 

DO  25  K = L , NPH I 
PH=PI2*(XIK|  + 1.) 

PFN=  PH*FN 
SN=S INC  « 5* PH) 

CRCK  ) = 4 . *SN * SN 
RI=PI2*A(K) 

Cl(K)  = .5*Rl*CR(K)*C0S(  PHN) 


C2CK  )=Rl*COS  CPH)+COS  (PHN  ) 
C3(K  )=R1*SIN< PH)*SIN(PHN) 
C4(K  )=P  HCnS(PFN) 

25  CONTINUE 
N2N=N2*N 
N4N=N2N  *2 


DO  62  J*  1 » N4N 

Y ( J ) =0  • 

Z C Jl *0 . " " 

62  CONTINUE 
U=(0.,  1.1 
DO  59  J= l » NG 
FJ=FN/RS< J) 

LI  =1 
L2=2 

IFCJ.LE.2)  L 1=2 
IFCJ.GT.N2)  L2  = l 
J8=J  + l 
J9=J8/2 
JT=2* J9* J-6 
J5=C J9-3 ) *N2-2 
SI  =1  . 

DO  60  I*1»J 
18=1*1 
19=18/2 
17=2*19+1-6 
J6= I 9+J5 
FI=FK/RSC I ) 

RP=R  SC  J )— R S C I ) 

ZP=ZSf JJ-ZSC  I» 

R2=RP*RP+ZP*ZP 

IF  ( I .NE.JI  GO  TO  41 

Sl  = . 5 

R2=.C625*D( J)*C(J) 

41  R3=R SC  1 1 *RS (J I 
G 1=0  . 

G2=0. 

G3=0. 

G4=0. 

G5=0. 

G6=0. 

CO  61  K= 1 » NPHI 
R4=R2+R3*CR (K ) 

R5=SCRTCR4) 

Zl=Sl/R5*CCCS(R5)-U*SI  NCR5H 
Y1=Z1*CUHI*R5)/R4 
G 1=C  1C  K»*Y1  + G1 
G2=C2(K1*Y14G2 
G3=C3CK1*Y1*G3 
GA=C  AC  K) *Z 1 *GA 
G5=C2C  K I *Z 1 +G5 
G6=C3CK»*Z1*G6 
6 1 CONTINUE 
G3  =U*G3 

YL=CRP*CVC  JI-ZP*SV(  jn*G2-RSCI»*CVC  J1*G1 

Y2=  C RS  C J )*SV C I)*CVCJ  I-RSCI)*SVC  J)*CVC  I )-ZP*SV(  IMSVC  JH*G3 
Y3  = Z P*G3 

Y4=  C RP  *CV( I >— ZF*SV( I ) ) *G2  + RS ( J)*CVC I)*G1 

Z1=U*CSV(IJ*SV (J)*G5*CVCI»*CVCJI*G4I 

G i=-U*G4 

Z2=-SV ( J )*G6 

G2=-FI*G4 

Z3=SVC  l ) *G6 

G3=F J*G4 

Z4=U*C G5-F I *G3 ) 


K2 =2 

IFII.LE.2)  K 1 = 2 
IFU  .GT.N2)  K 2-  1 
DO  3 i L =11 » L2 
LT=JT*L+L 
J7  = J6+L*N2 
DO  32  K*K 1 » K2 
KT= I T*K+K 
TT=T  (IT )*T ( KT  ) 

J 1= J 7*K 
J2=  J 1+N 
J3=J 1+N2N 
J4=J 3+N 

Y { J1  )=TT*Y1+Y  ( Jll 
Y(  J2)  = TT*Y2+YU2) 

Y{J3)=TT*Y3+Y(J3) 

Y<  J4 ) = TT*Y44YlJ4) 

ZIJ1  )=TT*ZI  + TP{LT)*TP(KT)*G1  + Zl  Jl) 

Z( J2  )=TT*Z2*TF(LT)*T (KT ) *G2  + Z C J 2 ) 

Z< J3)=TT*Z3+TF(KT)*T(LT)*G3*Z(J3) 

Z(J4  )=TT *Z4  + Z ( J4  I 
32  CCNTINUE 
3 l CONTINUE 
60  CONTINUE 
5«5  CONTINUE 
N2P=N2+1 
KD1  = l 
Jl=l 

00  II  J=1,N 
KD2=  KO 1*N 
KC3=KD1+N2N 
KD4=KD3+N 
R2=T  1J1  + 1I 
R3=T ( Jl+2) 

R4=T(JU3> 

J6=J5+1 

Rl=T(Jl)*T(Jl>*PD(  J5 )+R2*R2*  PD(  J6)  ♦R3*R3*PD(  J6+1 ) +R4*R4*P0  ( J6+2 ) 
G1=Y(K01HY(KD4» 

Y(KDI) =R 1+G 1 
Y(KD2)=C. 

Y t KD3 ) = 0 . 

Y ( KD4 ) =R l-Gl 
Z(KDi)=Z(KDl)*Z(KDll 
Z(KC2)=Z(KC2)-Z(K03) 

Z(KD3)=-Z(KD2) 

Z(K04I=Z(KD4)4Z(KD4) 

IF(J.EQ.l)  GC  TO  22 
KU1=KD  l-  1 
KU2=KD2-1 
KU3=KD3-1 
KU4=KD4- 1 
KLl=K01-N2 
KL2=KD2-N2 
KL3=KD3-N2 
KL4=  K04-N2 

Rl=T(J  l)*T( J1-2J*P0( J5)*R2*T( J1-1)*PD(  J6) 

G1=Y  (KUl»-Y(KL4) 

G2=Y(KU4)-Y(KU1 
Y ( KU l ) =R  i*G  i 
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Y(KU2I=Y(KL2)-Y(KL2) 
YIKU3l=Y(KU3)-Y(KL3) 
Y ( KU4 ) =P 1+G  2 

Y ( KL  1)  *R1-G2 
Y(KL  2 )=~Y(KU2  ) 
Y(KL3)=-Y(KU3  ) 

Y ( KL  4)  = R 1— G 1 
Z(KUl)  = Z<KUiMZ(KLll 
Z(KU2)=ZIKU2)-Z|KL3> 
ZIKU3)=Z(KU3)-Z(KL2) 
ZIKU4l=Z{KU4MZ(KL4) 
Z(KLI)*ZIKU1) 

ZIKL  2)=-Z<KL3 ) 

Z ( KL  3 ) =— Z ( KL2  1 
Z ( KL4 ) =Z I KU4 ) 

22  K0l=K01+N2P 
J 1=J  1+4 
J5=J5*2 
11  CONTINUE 

IFIN.LT.3)  RETURN 
J2-N2 

OC  13  1*3. N 
J2=J2*N2 
J1=I -2 
KL  t=  I 

00  14  J=l,Jl 
KU1= J2+ J 
KU2=KU1+N 
KU3= KU1+N2N 
KL4=  KU3+N 
KL2=KL 1+N 
KL3=KLl+N2N 
KL4=KL3+N 
Y(KL 1I=-Y(KU4 ) 

Y ( KL2 ) =-V I KU2 I 
Y ( KL  3) =-Y ( K03) 

Y ( KL4 ) = — Y ( KU 1 ) 
Z(KLI)=Z(KU1) 

Z ( KL  2 ) ~—H  KU3  ) 
Z(KL3)=-Z(KU2 I 
Z(KL4I=Z(KU4) 
KL1=KL1*N2 
14  CONTINUE 
13  CONTINUE 
RETURN 
END 


V.  THE  SUBROUTINE  PLANE 


A.  Description. 

The  subroutine  PLANE(NN,  N,  NT,  THR,  R)  stores  by  columns  in  R 
the  matrix  [R]  given  by 


ltQ 

RC* 

n 

n 

PQ 

n 

«*♦ 

n 

where  the  elements  of  the  column  vectors  on  the  right-hand  side  of  (12) 
are  given  by  (1-95) . The  input  variables  NN  and  N are  respectively  n 
of  (1-95)  and  the  maximum  value  of  i in  (1-95) . If  NT  = 1,  then  THR(l) 
is  0^  of  (1-95)  in  radians.  If  NT  > 1,  then  the  measurement  matrices 
[R]  of  (12)  for  6r  = THR(I) , I = 1,2, ...NT  are  stored  consecutively  in  R. 

The  variables  RS,  ZS,  SV,  CV,  and  T appearing  in  the  common  statement 
early  in  the  subroutine  PLANE  are  input  variables  calculated  by  calling 
any  one  of  the  subroutines  YMAT,  ZMAT,  or  YZ  beforehand.  The  calculated 
values  of  these  variables  depend  only  on  the  second,  fourth, and  fifth  arguments 
(NP,  RH,  and  ZH)  of  either  YMAT,  ZMAT,  or  YZ. 

Minimum  allocations  are  given  by 

COMPLEX  R(4*N*NT) 

DIMENSION  THR (NT) , BJ(M) 

COMMON  RS(NG),  ZS(NG),  SV(NG),  CV(NG) , T(4*N) 

where 

NG  = 2*N  + 2 

and  M is  the  largest  of  the  values  of  M calculated  by  PLANE.  The  sug- 
gested allocation  BJ(50)  will  work  if  the  maximum  circumference  of  the 
body  of  revolution  is  less  than  26  wavelengths. 

We  rewrite  (1-95)  as 
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i=2p-l 


i+2p-2 


where 

R. 


jkz  cos  9 


irj  (-2J  sin  0 cos  v+j  (.T  . . - J 1)cos  9 sin  v)e 
nr  n+1  n— 1 r 


- it  j (J  + J )cos  9 e 
J n+1  n-1  r 


jkz  cos  0 


ifj  (J  + J ,)sin  v e 
n+i  n-l 


jkz  cos  9 


. , Ikz  cos  0 

,n+l . T n r 

’3  (Jn+l  ' Jn-l,e 


In  (14),  p,z,  and  v are  to  be  evaluated  at  t = t^.  Equations  (13)  say 


i+2p-2Rl 

contributes 

to 

Rte 

np 

i+2p-2R2 

contributes 

to 

R*° 

np 

i+2p-2R3 

contributes 

to 

R^ 

np 

i+2p-2R4 

contributes 

to 

rH 

np 

ri  + 1, 

p - [—2 — ] - 1 


ri  + 1, 
p = [ ] 


Because  the  triangle  functions  do  not  overlap  at  the  ends  of  the  generating 
curve,  the  first  value  of  p must  be  discarded  if  i is  either  1 or  2 and  the 
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second  value  of  p must  be  discarded  if  i is  either  2*N4\L  or  2*N+2. 

In  DO  loop  12,  0^  of  (14)  is  THR(L) . DO  loop  13  evaluates  (14)  at 

t = tj.  The  logic  in  DO  loop  13  prior  to  statement  24  puts  the  Bessel 

functions  J , (x)  , J (x),  and  J . . (x)  in  BJ1,  BJ2,  and  BJ3  respectively 
n-1  n n+i 

where 


x = kpT  sin  0 
I r 

-8 

Choosing  M so  large  that  the  magnitude  of  JM_^(x)  is  around  10  , we 

start  with 


JM-l(x)  “ 0 
JM-2(x)  " 1 

and  use  the  recurrence  relation 

Jn-lW  ‘ T Jn<x>  - Jn+l^x)  <15) 

as  given  by  (9.1.27)  on  page  361  of  reference  [2]  to  calculate  Jn(x) 
for  n = M-3,  M-4,  ...  0 and  then  normalize  Jn(x)  according  to 

1 = JQ(x)  + 2J2(x)  2J4(x)  + 2J6(x)  + ...  (16) 

as  given  by  (9.1.46)  on  page  361  of  reference  [2].  Statement  24  and  the 
6 statements  following  it  put  R^,  R^,  R^»  and  R^  of  (14)  in  Rl,  R2,  R3 
and  R4  respectively.  DO  loop  20  adds  the  contributions  to  (13)  for  i = I 
and 


[2]  M.  Abramowitz  and  I.  A.  Stegun,  "Handbook  of  Mathematical  Functions," 
II.  S.  Govt.  Printing  Office,  Washington,  D.  C.  (Natl.  Bur.  Std.  U.  S. 
Applied  Math.  Ser.  55),  1964. 
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LISTING  CF  THE  SUBROUTINE  PLANE 


s 

! I 


f 

f 

f. 


I 


f 


SUBROUTINE  P L ANE <NN,  N, NT » THR « R ) 

COMPLEX  RI2401 ,U,Ul,U2 ,Rl,R2 »R3,R4 
,OIMENS  ION  TFR13I.BJI  501 

COMMON  RSI 42  I »ZS  142)  .SVI42I , CVI42)  tT  (80  ) 

N2=2*N 

NG=N2+2 

U=(0.,1.) 

U1=3.141593*U**NN 
JR=4*N*NT 
DO  22  J=  1* JR 
R(J)=0. 

22  CONTINUE 
J5=-2 

00  12  L = I t NT 
CS=CCS( THR { L ) ) 

SN=2.*SIMTFR(L)  ) 

DC  13  1*1, NG 
X=.25*RSU)*SN 
IFIX.LE..5E-7)  GO  TO  18 
M=2.8*X+ 13.-2. /X 
1FIX.LT..5)  M=10.8*AL0G10(XI 
IFIM.GE.INN+2) I GO  TC  19 

18  BJ1*C. 

8 J2=0. 

BJ3=C. 

IF (NN.EO.l ) ejl=l. 

IF(NN.EC.O)  BJ2=1. 

GO  TO  24 

19  BJt  M 1-0  . 

BJ(JM»=1. 

CC  14  J=3  * M 
JM=JM-1 

BJIJM)=JM/X*fiJ(JM+l I— B J ( JM+2 ) 

14  CONTINUE 

s=o. 

00  15  J=3 ♦ M , 2 
S=S+BJt  J) 

15  CONTINUE 
S=BJ(IH2.*S 
BJ2=BJ(NN*1)/S 
e J3=BJ(NN*2)/S 
BJl=-B J3 

IF(NN.GT.O)  PJ1=BJINN)/S 
24  ARG=ZSCI)*C$ 

U2=U  1*1  CCS ( ARG ) +0*Sl N( ARG)) 

P4= ( EJ3-EJ l )*U*U2 
R2  = ( BJ3+BJl)*U2 

Rl=-PJ2*CV( I >*$N*U2+CS*SVU )*R4 

R3=SV( I ) *R2 

R2=-CS*R2 

19- ( 1 + 1 )/2 

I T=2  * 1 9*  1-6 

J7=I 5+ J5 

Kl  = l 

K2=  2 

IFII.LE.2)  K 1 = 2 
IFU.GT.N2)  K2=l 
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00  20  K=K  1 » K 2 
TV=T(lTtK+K) 

J 1= J 7+K 
J2=J1+N 
J3=  J 2+N 
J4= J3+N 

R(  J1  )=TT*Rl*R( J1 ) 
R ( J2  J=TT*R2  + R { J2) 
R(J3J=TT*R3+R(J3) 
R(J4)  = TT*R4'»-R(J4) 
20  CONTINUE 
13  CONTINUE 
J5=J<t-2 
12  CONTINUE 
RETURN 
END 
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VI.  THE  SUBROUTINES  DECOMP  AND  SOLVE 


A.  Description; 

The  subroutines  DECOMP (N, IPS, UL)  and  SOLVE (N, IPS, UL,B,X)  solve  a 
system  of  N linear  equations  in  N unknowns.  These  subroutines  will  be 
used  in  Section  VII  to  solve  the  H-field  matrix  equations  (1-17),  the 
E-field  matrix  equations  (1-57)  and  the  combined  field  matrix  equations 
(1-88) . The.  input  to  DF.COMP  consists  of  N and  the  N by  N matrix  of 
coefficients  on  the  left-hand  side  of  the  matrix  equation  stored  by 
columns  in  UL.  The  output  from  DECOMP  is  IPS  and  UL.  This  output  is 
fed  into  SOLVE.  The  rest  of  the  input  to  SOLVE  consists  of  N and  the 
column  of  coefficients  on  the  right-hand  side  of  the  matrix  equation 
stored  in  B.  SOLVE  puts  the  solution  to  the  matrix  equation  in  X. 

Minimum  allocations  are  given  by 

COMPLEX  UL(N*N) 

DIMENSION  SCL(N) , IPS(N) 

in  DECOMP  and  by 

COMPLEX  UL(N*N),  B(N) , X(N) 

DIMENSION  IPS(N) 

in  SOLVE. 

A description  of  DF.COMP  and  SOLVE  is  on  pages  46-49  of  reference  [3]. 


[3]  J.  R.  Mautz  and  R.  F.  Harrington,  "Transmission  from  a Rectangular 
Waveguide  into  Half  Space  Through  a Rectangular  Aperture,"  Interim 
Technical  Report  No.  12,  Rome  Air  Development  Center,  Griff iss  Air 
Force  Base,  New  York,  August  1976. 
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listing  cf  tee  subroutine  oecomp 

SUBROUTINE  CECCMPIN. IPS.UL) 

COMPLEX  UL ( 160C I tP IVOT *£M 
DIMENSION  SCL(40),IPS(40) 

00  5 1=1. M 
IPS!  l)  = I 
RN=0  . 

Jl  = I 

ULM*ABS(REAL (UL ( Jl H J+ABSIAIHAGIUU  Jl)  > > 

Jl=JI+N 

IF(RN-ULM)  1.2.2 

1 RN=ULM 

2 CONTINUE 

SCL ( I )= 1 /RN 
5 CONTINUE 
NM 1=  N— 1 
K2=0 

DO  17  K = t » NN1 
BIG=0. 

DO  11  I=K,N 
IP=IPSU) 

SUEMABSIREAUULIIPK)  ) | ♦ASS  (AIMAGCUU  IPKIH  l*SCL  ( IP » 

IF ( S 12  E-BIG ) 11.11*10 

10  B IG=S I ZE 
IPV  = I 

11  CONTINUE 
IF(IPV-K)  IA.15.14 

14  J=IPS(K) 

IPS(K)=  IPS ( IPV) 

IP S ( IPV)=J 

15  KPP=IPS(K)+K2 
P!VOT=UUKPF> 

KP1=K+ 1 

DO  16  I=KP1.N 
KP=KPP 

IP=IPS(I1+K2 
EM=-UL ( IPl/F IVOT 
L 8 IJL(IP)=-EM 
no  16  J=KPi,N 
l P= I P+N 
KP=K  P+ N 

UL(  1 P)=UL( IP  ) + FM*UL(KP  ) 

16  CONTINUE 
K2=K2*N 

17  CGNTINUE 
RETURN 
ENO 

LISTING  OF  THE  SUBROUTINE  SOLVE 

SUBROUTINE  SCLVEIN. I PS.UL. B» X) 

COMPLEX  UK  16CC ) *8(40)  »X(40)  .SUM 
CIMENSION  I P S( 40) 

NPKN+l 
1P=IPS( U 
Xtl)=B( IP) 

DO  2 i = 2*N 
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IP=IPS< I ) 

I PB=  IP 
IP1= !-l 
SUM=C. 

DO  1 J = 1 * I M 1 
SUH=$tjMUL( IPI*Xl  J) 

1 [ P» I P+N 

2 X(  I)=BUPB)-SUP 
K2=N*(  N-  l) 

IP=  l PS  I N ) +K  2 
X( N> =X(N)/UL<  IP) 

CO  A IBACK=2,N 

I"*NPI-  1 8ACK 

K2=K2-N 

IPl  = IPS 1 1 ) +K  2 

IP1  = U1 

SUM=C. 

I P-  I P I 

DO  3 J=IPltN 
IP= I P+N 

3 SUP=SUM*UL( IP)*X(J) 

A >n  ) *1  xm-SLM  ZULU  PI  ) 
BETUPN 
END 
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VII . THE  MAIN  PROGRAM 


A.  Description: 

The  main  program  uses  all  the  subroutines  listed  in  Sections  II-VI 
to  evaluate  the  surface  currents  (1-109)  at  one  specified  value  of  <j>  for 
both  0 and  <f>  transmitter  polarizations  and  to  evaluate  the  scattering  crocs 
section  (1-111)  at  one  specified  r1  • 3 of  <fr  for  the  four  different  combi- 
nations of  transmitter  and  receiver  polarizations.  The  main  program  is  not 
general  enough  to  meet  every  user's  needs.  Its  main  purpose  is  to  give  an 
example  of  how  these  subroutines  can  be  used  and  to  provide  the  use:  with 
sample  output  from  them  so  that  he  can  be  sure  that  they  are  working 
properly. 

Input  data  is  read  from  punched  cards  according  to 

RE AD (1,51)  NM,  NP,  NPHI 
51  FORMAT  (313) 

READ(1,50)  BK,  TT,  P,  TR,  PR,  ALP 
50  FORMAT  (5E14.7) 

READ(1,46)(RH(I),  1=1,  NP) 

READ(1,46)(ZH(I),  1*1,  NP) 

46  FORMAT  (10F8.4) 

READ(1,50)(X(K),  K = 1,  NPHI) 

RF,AD(1,50)  (A(K),  K » 1,  NPHI) 

The  above  input  variables  are  defined  by 


NM  » one  plus  the  maximum 

value  of  n in  (1-109) 

NP  = P 

NPHI  = N, 

9 

BK  = k 
TT  = 0 

t 

P = 4> 


tr  = e 

r 

and  (1-110)  PR  = * 

ALP  » a 
RH(i)  - kp~ 
ZH(i)  = kz" 

X(k)  * x^ 
A(k)  = A^ 


where  all  variables  on  the  right-hand  sides  of  the  above  equations  appear 

In 

in  reference  [1].  According  to  reference  [1],  n denotes  eJl  v dependence, 
(p i , z±),  i=l,2,...P  are  coordinates  on  the  generating  curve,  k is  the 
propagation  constant,  0 is  the  colatitude  of  the  transmitter,  and  <j>  is 
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the  longitude  at  which  the  surface  current  is  evaluated.  The  angles  6^  and 
$ are  the  colatitude  and  longitude  of  the  receiver.  The  real  constant  a is 
the  relative  weight  of  the  E~field  integral  equation  in  the  combined  field 
formulation.  Finally,  x^  and  are  the  abscissas  and  weights  for  the 
point  Gaussian  quadrature  integration  in  $.  All  angles  read  into  the  main 
program  are  in  degrees. 

Minimum  allocations  are  given  by 

COMPLEX  TJ (6*N) , FJ(6*N),  RT(4*N),  RR(4*N) , 

Y (4*N*N) , Z(4*N*N),  B(2*N),  C(2*N) 

DIMENSION  RH(NP),  ZH(NP) , X(NPHI) , A(NPHI), 

THT(l) , THR(l) , R2(N),  IPS(2*N) 

where 

N *=  (NP  - 3)/2 

The  surface  current  printed  by  the  main  program  is  not  (1-109)  per  se, 
but  (1-109)  divided  by  the  propagation  constant  k.  The  heading  above  the 
surface  current  in  the  printed  output  denotes  real  and  imaginary  parts  and 
t and  <J>  components.  As  further  identification,  NHEC=1  denotes  the  H-field 
solution,  NHEC=2  denotes  the  E-field  solution,  NREC=3  denotes  the  combined 
field  solution,  KT  = 1 denotes  a 0 polarized  transmitter,  and  KT  ® 2 denotes 
a <f>  polarized  transmitter.  The  scattering  cross  section  per  square  wave- 
length (1-111)  is  labeled  SIGMA/ (LAMBDA) **2  and  further  identified  by  the 
parameters  NHEC,  KT,  and  KR  where  KR  « 1 denotes  a 0 polarized  receiver  and 
KR  = 2 denotes  a $ polarized  receiver. 

The  index  K of  DO  loop  41  obtains  the  n*=r-l  term  in  (1-109)  and  (1-110). 

In  DO  loop  58,  NHEC=1  obtains  the  H-field  solutior , NHEC=2  obtains  the  E-field 
solution,  and  NHEC=3  obtains  the  combined  field  solution.  In  DO  loop  27, 

KT®!  denotes  a 0 polarized  transmitter  and  KT*=2  denotes  a $ polarised  transmitter. 
In  DO  loop  16,  KR®1  denotes  a 9 polarized  receiver  and  KR®2  denotes  a 4>  polarized 
receiver.  The  user  can  omit  computations  which  are  of  no  interest  to  him  by 
changing  the  statements  which  introduce  DO  loops  41,  58,  27,  and  16  so  as  to 
restrict  the  indi..'  <"f  these  DO  loops.  The  logic  in  DO  loop  58  before  state- 
ment 59  puts  the  square  matrix  on  the  left-hand  side  of  either  the  H-field 
matrix  equation  (1-17),  the  E-field  matrix  equation  (1-57),  or  the  combined 
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field  matrix  equation  (1-88)  in  Y.  The  logic  in  DO  loop  27  before  state- 
ment 53  uses  (1-99),  (1-100)  and  (1-104)  to  put  the  excitation  column 
vector  on  the  right-hand  side  of  either  the  H-field,  the  E-field,  or  the 
combined  field  matrix  equation  in  B.  Statement  53  puts  the  solution  of 
the  matrix  equation  in  C.  If  KT=1,  the  logic  in  DO  loop  13  adds  the 
n“K-l  term  of  the  t and  4>  components  of  the  numerator  of  the  first  of 
equations  (1-109)  to  TJ  and  PJ  respectively.  DO  loop  15  adds  the  n=K-l 
term  of  the  t and  <j>  components  of  the  numerator  of  the  second  of  equations 
(3-109)  to  TJ  and  PJ  respectively.  DO  loop  16  adds  the  n=K-l  term  of 


4irr 


r 

-jn 


<J) 


to  F.(KE)  where  E^(J)  is  given  by  (1-110).  Here,  p is  either  0 or  and 
q is  either  0 or  4>. 

Nested  DO  loops  28,  29,  and  30  print  out  both  (1-109)  divided  by  the 
propagation  constant  k and  (1-111).  The  indices  NHEC,  RT,  and  KR  of  these 
DO  loops  have  the  same  meaning  as  in  nested  DO  loops  58,  27,  and  16. 

The  sample  input  and  output  data  is  for  oblique  incidence  on  a 
conducting  sphere  of  radius  0.2  wavelengths.  Only  the  n=0  and  n=l  terms 
in  (1-109)  and  (1-110)  are  considered.  A few  more  terms  may  be  needed  to 
obtain  reasonable  accuracy. 
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B.  LISTING  OF  THE  MAIN  PROGRAM 

//PGM  JOB  (XXXXfX>XX*2»2)t' MAUTZ, JOE*  , REGI0N-2OOK 
//  EXEC  WATFIV 
//GO.SYSIN  OD  * 

I JOB  MAUTZ, TIME=2,PAGES=40 

C THIS  PROGRAM  RAN  FOR  45  SECONDS  ON  THE  IBM  370/155. 

C SLBRCUTINE  S YMAT,  ZM AT , YZ,  PLANE,  DECQMP,  AND  SCLVE  ARE  CALLEC. 

COMPLEX  U,TJ(I20),PJ(I20),E(  12 ) , SN , SNR  ,RT(  240)  ,RR(  240  , Y ( 1 fcOC) 
CGMPLEX  2(1600) ,B( 40 ) ,C (40 ) , U l,CCNJG 

DIMENSION  RH(42) ,ZH(43) ,X(20) ,A( 20* ,THT(3) ,THR(3) ,R2(20> , I PS (40) 

R E AD (I ,5 1 ) NM,NP,NPHI 
51  FORMAT (313) 

RE  AO ( I, 50  ) BK, TT,P,TR,PR, ALP 
50  FORMAT ( 5E14 .7  ) 

READ ( 1,46) (RH( I ) ,1=1 ,NP» 

READ  Cl  ,46  M ZM  I),I=l,NP> 

46  F ORM AT ( 1 OF 8.4 1 

RE AD ( l » 50 ) ( X ( K I ,K= 1 , NPHI ) 

R6A0(l,50)(A(K),K-l,NPHII 
WRITE(3,49)  N M ,NP , NPHI 
4D  FORM  AT ( • NM  NP  NPHI •/ IX, 21 3, 141 
WB I TF ( 3 ,48 ) BK,TT,P,TR,PR,ALP 

4 8 FORM  AT l 7X, *EK« ,I2X, 'TT • , 13X,§ P» , 12X, »TR» , 12X , • PR  * />  X ,5E 14 . 7/7X , • AL 
1P*/IX,E14.7) 

WRITE(3,45) ( RH ( I ) , 1= 1 , NP) 

45  FORM  AT  C • RH»/ ( IX, 10F8.4)  ) 

WRITE(3,44) (ZH( II,I=1,NP) 

44  FORMAT! • ZH • / ( 1 X , IGF  8 .4 )) 

WRITE(3,47) (X(KJ,K=l,NPHl| 

47  FORMAT ( * X'/(1X,5FI4.7») 

WRITE! 3, 43) (A(K) ,K=l,NPHl ) 

43  FORMAT ( • A* / ( IX, 5E14 .7 ) * 

N2=NP-3 

N=N2/2 

N3=N2*N 

N4=N2*N2 

N6=6*N 

U=(0.,I.) 

P 1=3 .141593 
P4=. 062 5/P  I** 3 
P 8=P  I / 1 8C. 

THT ( l )=TT*P8 
P = P*P8 

THR(  l)  = TR*Pfi 
PR=PP*P8 
DC  42  J* 1 , NP 
RH ( J )=BMRH(J  ) 

ZH( J)=BK*ZH( J) 

42  CONTINUE 
CO  17  J*1,N 
R2(J)-I./RH(2*J+I) 

17  CONTINUE 

00  54  J=1,N6 
T J ( J ) =0. 

PJ(J»=0. 

54  CONTINUE 

DO  55  J=  I , I 2 
F ( J)  *0  . 

55  CONTINUE 
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WRITE(3,9) 

9 FORMAT ( * OSA^PLE  OUTPUT  FROM  SUBROUTINES*) 
DO  41  K=1,NM 
NN=K-1 
PN=NN*P 
C$=CDS (PN) 

SN=2  •* SINI  PM  *L 

PN=NN*PR 

CSR=COS!FN) 

SNR=2.*SIN<  PM*U 
IF(NN.EC.O)  GO  TO  56 
C S=2  .♦CS 
CSR=2.*CSR 

56  LANE=0 

DO  58  NHEC  = 1,3 

GO  TO  ! 6 1,  6 2,  £ 2 ) » NHEC 

61  CALL  YMAT!NN,NF,NPHI,RH,ZH,X,A,Y) 

WRITE! 3t8)  Y(1),Y<2) 

8 FORM  AT !1X»4E14»7) 

GC  TC  59 

62  CALL  ZMAT! NN,NP,NPHI ,RH, ZH, X,A,Y> 

V,  R I T E ( 3 1 8 ) Y11),Y12) 

GO  TO  59 

63  CALL  YZINN* NP»NPHItRHtZH»X.A,YfZ) 
WRITE(3f8)  Y{ 1 ) ,Y(2)  ,ZI1) tZ(2) 

00  66  J=1,N4 
Y(J)=Y(J)+ALP*2lJ) 

66  CONTINUE 

59  CALL  DEC0MPIN2,IPS,Y) 

WRITE!3,8)  Ytl),Y(2) 

IF(LANE.NE.C)  GO  TO  57 
LANE=1 

CALL  PLANE! NN»N»i.THT,RTI 
CALL  PLANE!NN,N,ItTHRfRR) 

WRITE  !3 ,8)  RT(I),RTC2) ,RRlll,RRI2> 

57  DO  27  KT=l,2 
L=2*  INHEC-I )*KT 

GO  TC  131, 32, 33, 34, 35,36), L 

31  DO  21  J= 1 » N 

B ! J ) =-RT 1 J +N  3 ) 

B! J+N) =— RT { J+N2 ) 

21  CONTINUE 
GO  TO  53 

32  DO  22  J= 1 » N 
JN= J +N 

B!  J) =- R T 1 J N ) 

B1JN  )=-RT( J ) 

22  CONTINUE 
GO  TC  53 

33  DO  23  J=  1, N 
B( J)=RT( J) 

JN= J+N 

B ( JN )=-RT ( JN  ) 

23  CONTINUE 
GO  TO  53 

34  00  24  J=1,N 

B!  J)  =— RT ( J+N2  ) 

B( J+N )=RT 1 J*N3  ) 

24  CONTINUE 
GO  TO  53 
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35  00  25  J*  l , N 

B(  J)  =-RT( J+N3)+ALP*RT(  J) 

JN=J*N 

B< JN)*-RT« J+N2  )-ALP*RT 1JN) 

2 5 CONTINUE 
GO  TO  53 

36  DO  26  J = i » N 
JN=J*N 

B < J)  =-RT (JN)-ALP*RT( J+N2 ) 

B ( JN )=-RT( J ) *ALP*R  T I J+N3) 

26  CONTINUE 

53  CALL  $CLVE!N2,IFS,Y,B,C) 

WRITE! 3»  8)  C(1),C(2) 

Jl  = ! L-l ) *N 
GO  TO  I 11*12) » KT 

11  DO  13  J»1,N 
J2=J+J1 

TJ( J2)=TJ< J2)+C(J)*CS 
PJ(J2)=PJ<  J2MC(J+N)*SN 

13  CONTINUE 
GO  TO  14 

12  DO  15  J = 1 » N 
J2=J+Jl 

TJ(J2)«TJ(J2MC(J)*SN 
PJI J2) =PJ( J2)+C!J+N)*CS 

15  CONTINUE 

14  DO  16  KR=1,2 
L=( KT- 1 ) *2*  KR 
KE=4*(NHEC-1)+L 
U1  = 0. 

GO  TO  (71, 72,73,74), L 

71  DO  75  J=l,N2 
Ui=U  1+RR  ( J)  *C  l J) 

75  CONT  INUE 
E(KE)=E(KE)+U1*CSR 
GO  TO  16 

72  DO  76  JM,N2 
Ul=Ui*-RR(J  + N2)*C(J) 

76  CONTINUE 

E(KE  )=E(KE)*Ui*SNR 
GO  TO  16 

73  DO  77  J=  1, N 2 

U1  = U1  + RR  !J)*C ! J) 

77  CONTINUE 

EIRE  )=E(*E  )*Ul*SNR 
GO  TO  16 

74  DO  78  J*l,N2 
Ul=Ul«-RRIJ  + N2)*C(J) 

78  CONTINUE 
E(KE)=E!KEHU1<CSR 

16  CONTINUE 

27  CONTINUE 
58  CONTINUE 
41  CONTINUE 

DO  28  NF  EC= 1,3 
DO  29  K T=1  , 2 
WR ITE! 3, 18)  NFEC  ,K  T 

18  FORMAT ( *0NFEC=*,I3,' , KT=*,I3I 
WRI TE! 3,  19) 

19  FORMAT!*  REAL  JT  IMAG  JT  REAL  JP 


IMAG  JP* ) 
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Jl  = N*t2*f\HFOKT-3) 

DO  37  J=1,N 
J2=U+ J1 

TJ( J2)=TJ( J2)*R2< J) 
PJ1J2)=PJ(J2)*P2(J) 

WR I T F.  ( 3 » 38  ) TJ(J2),PJ{ 
38  FORM  AT ( IX, 4E 11.4) 

37  CONT  INUF 

DO  3 C KR  = 1 1 2 
J1=4*NHEC+2*KT+KR— 6 
S!G=P4*E(Jl )*CCNJG ( E ( J 
WRITE( 3 , 10)  NEEC  , KT, KR 
10  FORMAT  (•  NHEC=  * » I3»  • , 
11.4) 

30  CONT  1NUE 
29  CONTINUE 
28  CONTINUE 
STOP 
END 

! AT  A 

2 21  20 

: . 1 256637F+  0 1 C. 3 C C C CCOE+  02 
: .2  000000  E+00 

0.0000  C . 1 564  0.3090  0. 

1.0000  0.9877  0.9511  0. 


1>> 

t SIG 

KT= ' ♦ I 3* • f KR=  • , 


13, •»  SIGMA/(LAMBDA)**2=‘  ,E1 


C. 3CCCCC0E+02  0.4500000E+02  0. 60000008+02  0 .20000C0 E+02 


0.3090 

0.9511 

-0.9511 

0.3090 


0.0000  C . 1 564  0.3090  0. 

1.0000  0.9877  0.9511  0. 

C.0000 

-1.C3C0  - C.  9877  -0.9511  -0. 

0. 0000  0.1564  0.3090  0. 

1.  COOO 

:. 9931286  E+OO-C. 963  9 719E+C0 
) .6 360537 E+00-0.51C8670E+00 
).7652652E-01  0. 2277859E+00 
! .7 463319E+00  0.839  l 17GE+G0 
) . I 76 1401 E -0 1 0.4060143E-01 
:.l  181945E+00  0.  1316886E+C0 
) .1527534E  + 00  0 . 149  1730E  + 00 
5 .1 019301  E + 00  0 . 832  7674E-01 


! .7 463319E+00 
i . I 76 1401 E -0 1 
.1 181945E+00 


4540  0.5878  0. 

8910  C . 8090  0. 

8910  -0.8090  -0. 

4540  0.5878  0. 

-C. 91223448+00-0 

-Q.3737061E+00— 0 
0.3737061 E+OQ  0 
G.9122344E+00  0 
0.6267205E-01  0 
0.142G961E+00  0 
0.142C961E+00  0 
0.626  7205E-0 1 0 


7071  0.8090 

7071  0.5878 

7071  -0.5878 
7071  0.8090 

•8391170E+00- 
• 2277859E+  CO- 
•5108670E+00 

• 96397 19E+  00 
•8327674E-01 
• 149 1730E+0G 
. 13 16866E+  CO 

• 40601 43E-0 1 


C . 85 10 
0.4540 


C .95 1 1 
0.3090 


-C .4540  -C.3090 
0.8910  0.9511 

0 .74633196+00 
0 • 76  526  52E-C1 
0.6360537E+00 
0.9931286E+00 
0.10  19301E+00 
0 • 1527534  E+00 
0.1 181945E+00 
0 . 17614C IE— 0 1 


FF  INTEO  OUTPUT 
N*  NP  NPH I 
2 21  20 
BK 

C.l 256637E+01 
ALP 

C .2 JJ JCCOE+OO 
RE 

0. 0000  0.1564 

1.  0300  0.9877 

C.0000 

7H 

- i.  0000  -C.98 
0.0000  0.15 

1. 0000 


C. 9 877 
0.1564 


TT  P TR  PR 

3CCC0GGF+02  0 .45QC000E+02  0 .60000008*02  0 . 20000 C0E+02 


0.3090  0.4540  C.5878  0.7071  0.8090  C.8910  0.9511 

0.9511  0.8910  C. 8090  0.7071  0.5878  C.4540  C.3090 


■C.5511  -0.8910  -0.8090  -0.7071  -0.5878  -0.4540  -C.3090 
0.3090  0.4540  0.5878  0.  7071  0.  8090  C.8910  C.95U 


■C  .993  1286E+00-C 
•C.6360537F+00-C 
C .7652652E-CI  C 


.9fc2S719E+0Q- 0.912  2344E+00— 0. 8391170E+  00- C .7463319E+00 
,5108670 E + 00-0. 373 7061 E+00-0 . 2277859 E + 00-0 .765265 2 E-Cl 
,22776  59E * 00  0.3737061 E+OQ  0 . 5 1 0 86 70E  + 00  0 .6360537E+00 
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C.7463319E+00  0.8391170E+00  0.9122344E+00 
A 

0.1761401E-01  0 .40601436-01  0 .62672C8 E-01 
0.1  181945E+00  0.  1316886E+C0  0. 142 C961 E+00 
0 .15275346+00  0. 149  1730E+0Q  0. 14209616+00 
C . 1 019301 E+00  0.8327675E-01  0 .6267208 E-0 l 


0.96397 19E+00 

0.8327675E-CI 
0. 14917306+00 
0. 131 68  86E  + CO 
0 .4060143E-Q  1 


SAMPLE  CUT  PUT  FROM  SLPROUTIN 
0.2  0969146+01  0 . 1 786422 E-Ol 
C.2096914E+01  C . 1 78*422E-C1 
-C.4 120C92E-01  0 .9  162729E  + 00 
- G .42436536+  00  G. 5C08702E+00 
-0.1092077E+00-0 .10348446-01 
C.OOOOOOOE+OO  O.CCCCOOOE+OO 
C.3091431E— 01— 0. 1C  7 5599E+  02 
C .309 1431  E-Ql— 0 . 10  755996  + 02 
-0.1 23 0 594 E+ 00-0 . 1 179976E-01 
O.OOOOOOOE+OO  O.CCOCOOOE+OO 
C.2096915E+01  0 . 1786922E-01 
C.3091431E-C1— 0. IC75599E+02 
0.2 103098 E+0 1-0.2 133332E+01 
-0.1 184121 E+00-0 .4 237223E-02 
O.OOOOOOOE  + OO  C.CCOCOOOE+OO 
0.2 437500 E+01 -0.10570376-01 
C.2437500E+01-0. 10570376-01 
C.4 8071186*00— 0.87 3 1C53E+00 
C. 3939573 E+00-0. 4160154E+03 
C.8035219E-01  0 . 2 195775E+00 
0.1 91 8676E+00-0. 13720976+00 
C • 140 8945 E+  00-0. 5 3 85 362 E+ 01 
C. 140 8945 E+00-0. 5 385362E+01 
C • 8862776E-01  0 . 2205508E+03 
C.1875597F+00-0. 1473804E+00 
0, 2437500 E+0 1-0. 1057036E-01 
C, 1 40 8945 E+OC— 0.5389362 6+ Ql 
C .2 46 56 78 E+0 1-0 . IC8844  2E  + 01 
C • 7 738441 E-Oi  0 .2207815E  + 00 
C.1917793E+00-0. 1352912E+C0 


ES 

0.28 27 5 16 E+00  0 
C. 134961 7E+00  0 
-0.2257579E+00  0 
-C. 84 16 800 E+00  0 
-0.52914926+00-0 
O.OOOOOOOE+OO  0 
C.5109501E— 01  0 
-0.65065626-01  0 
-0. 51 6 2 89 8 E+00— 0 
0. OOQCOOOE+OO  0 
0.28275 14E+00  0 
0.5109501E— 01  0 
C.2778168E— 01  0 
-0 .52 16 5 29 E+00— 0 
O.COOCOOOE+OO  0 
0 .51409 11 E+00- 0 
0 .21094386+00-0 
0.420 1249E+00-C 
0. 1906640E+00-0 
0.7623661E— 02  0 
0 .332 95 19 E+00-0 
0* 1184365E+00  0 
-0.305 2 041 E+00  0 
0.5910277 E-02  0 
0. 3282 781 E+00— 0 
0.5 1409 10 E+00-0 
0.1184365E+00  0 
0. 1363794E+00  0 
0.505 1553E-02  0 
0. 33 1759 6E+ 00-0 


•3195152E-01 
.1408762 E-Oi 
•7185605E+C0 
.78531276+00 
.10062126+00 
.00000006+00 
.69999 79E+00 
•49373846-C2 

• 1013399E+00 
. OCOOOOOE+OO 
•3195150E- 01 
.6999979E+00 
. 1099419E+CC 
.1052166 E+00 
.OCOOOOOE+OO 
•2176931E-01 
•8016225E-02 

• 7 168542E+  CO 
.54760636+00 
.35301516+00 
•9042 1 14E- 0 1 
• 1649075E+01 

• 2995496E—  01 
•3503525E+00 
. 8860910E— 01 
.2176930E-C1 

• 1649075  E+01 

• 1851363E+  00 
.35204486+00 

• 8850223E— 01 


NF  EC=  1,  KT=  1 

REAL  JT  IMAG  JT  REAL  JP  IMAG  JP 
0.1 140E-01  0 .774 5 E+00-0 .30956+00-0 .7744 E+00 
-C.7018F+00  0.5357E+CO-0.2263E+00-C.5920E+00 
-0.1  187E+01  0.1405E  + CC-0.2110E+00-0.3438E  + 00 
-C. 1 440E+01-0 .353  8E  + C0-0 .3 40 8 E+00-0 .1353E+00 
-C. 1 41 0E+01-0.8C99E  + C0-0.5889E+00-C. 10 086+00 
-0.1 199E+0 1-0.1 133E+01-0.8275E+00-C.2960E+00 
-C.9433E+00-0.  13C8E  + C 1-0.9  31  9E  + 00-C.5'>89E+00 
-0 .7  49 8 E+0 0-0 .1383E  + C 1—0.8  846E  + 00- O.iOlOE  + Oi 
-C.7 116E+00-0. 13 72 E+C l -0.7766E+00-0 .12666+01 
NF  EC  = 1.  KT=  1,  KR=  1,  SI GMA/ ( L AMBD A)**2=  0.3550E+00 
NF  EC-  1*  KT  = 1,  KR=  2,  S I GMA/ (L AM  BOA) **2=  0.3835E-01 


NF  EC  = 1«  KT=  2 

PEAL  JT  IMAG  JT 
C.4997E+00  0.6988E+C0 
C .1 731 E+00  0 .6375  E + CC 
-C.2673E+00  0. 480 1E+00 
-0.7001E+00  0.2035E+CC 
-0.1004E+01-O.1645E+C0 


REAL  JP  IMAG  JP 
0 .51006+00  C .52746+00 
0.58126+00  0 • 1925E+00 
0.84146+00  0.6892E-01 
0 . 1 144E  + 01  0. 2235E  + 00 
0 ♦ 1 288E+01  0.6264E+00 


.99312866+00 

• 1C193C1E+C0 
.15275346+00 
. U81945E+00 
•17614C1E-01 
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-0.1116E+01-0.5484E+CC  0.1190E+01  0.1076E+01 
-0.1058E+01-0.8667E+CO  0.9518E+00  0.1381E+01 
-C.9129E+00-0. IC75E+C1  0.7419E+00  C.1496E+01 
-0.7720E+00-0.1177E+C1  0.6952E+Q0  C.1346E+01 


NH6C  = 

1. 

KT  = 

2, 

II 

a 

It  SIGMA/ ( L AM80A )**2= 

NE  FC  = 

I* 

KT  = 

2, 

KR  = 

2 « SIGMA/(LAMBI3A)**2= 

NE  EC  = 
REAL 

2, 

JT 

KT=  1 
IMAG 

JT 

REAL  JP  IMAG  JP 

C .5  869E-02  0.7729E+00-0.3745E+00-C.8135E+00 
- C .6  876E  + 0C  0.5336E+C0-Q.  1304E+00-C.5761E+00 
-0.1 194E+01  O.1322E+C0-O .2 399E+00- C.35 14P+0G 
-C.1434E+01-0.3588E+CO-3.3084E+00-U.,  j24E-0l 
-C.1408E  + 01-0.8132E  + CO-0.604  7E  + 00-Q.8353E-C1 
-C .1  191E+0 1-0.1136E+CI-0.8651E+00-C.26G6E+00 
-C.9  372E+OC--0.  13  12E+  C 1-0. 9 41 1 E+00-C  .67  77E+00 
-0.7  38  2E+00-0.1390E  + 01-0.9  73  2E  + 00-0.10UE  + 01 
-C.6682E+00-0.1430E+01-0 .7495  E 400-0 • 1341 E *01 
NEEC  = 2,  KT  = 1,  KR  = 1,  SI GMA/ C LAMBDA! **2- 
NE  EC  = 2 ♦ KT  = 1,  KR=  2,  S I GMA/ (L AM  80 A ) **2= 


0.2916E-01 

0.3866E+00 


0.3546E+0G 

0.3856E-01 


NE  EC=  2,  KT=  2 
REAL  JT  I RAG  JT 

0.5368E-+00  0 .6831E  + C0 
0.1696E+00  0.6285E+00 
-C.2848E+0C  0.4730E+CO 
-0  .7 1 16E+0 0 0 • 19  75E  + 00 
-0.1017E+01-0.1681E+00 
-0.1122F+01— 0.5508E+C0 
-0.106  5E+01— 0.8633E  + C0 
-C.9066E+0C- J. 1079E+01 
-0.7407E+0 0-0.1 18 1E+C1 
NH  EC  = 2,  KT=  2,  KR  = 
NHEC  = 2 1 KT,  2,  KP  = 


REAL  JP  IMAG  JP 
0.5764E+00  C.5256E+00 
0.4836E+00  C • L829E+00 
0.8914E+0Q  C.3695E-01 
0.1144E+01  C. L921E+00 
0.1335E+01  0.6229E+00 
0. 1233E+01  C.1095E+01 
0.948 IE+00  C • 1433E+01 
0. 8055E+00  0. 1541E+01 

0 . 6, 582E  + 00  C. 1416E+01 

1,  S I GM A/ (L  AMBOA )**2= 

2,  SIGMA/<LAMBDA)**2= 


0.2954E-01 

0.3880E+00 


NE  EC  = 3,  KT=  1 

REAL  JT  IMAG  JT  REAL  JP  IMAG  JP 
-0.231 1 E-Oi  0.7932E+00-0.3317E+00-C.77 18E+00 
- C .6 966E  + 0 0 0.53 16E+CG— 0.2239E+00-C.5899E+00 
-0.1  19  5E+01  0.1346E+C 0-0.211 5E+00- 0.341 7E +00 
-G.1438E+01-0.361 7E+C0- 0.340 IE+00— 0 .1333E+00 
-C .1 406E+01-G.8171E+CO— 0.5890E+00-C.9960E-01 
-0.1 191E+01-0.1138E+01-0.8276E+00-C.2952E+00 
-C.9329E+0C-0. 13 1 1F+C1-C.9335E+00-C.6489E+00 
-0.7  433E+0  0-0. 13  78  E + C 1-0.8  864E+00-0. 10  07E  + 01 
-C.7Q68E+00-0.1370E+CI-0. 7 994  E +00— 0 • 1275E+Q 1 
NE  EC  = 3,  KT=  It  KR  = 1,  $IGMA/UAMBDA)**2=  0.3546E+00 
NE  EC  = 3,  KT  = I,  KR=  2,  SI GMA/ (L AM80A )**2=  0.3821E-01 


NE  EC  = 3,  KT=  2 

REAL  JT  IMAG  JT 
C .4927E+00  0 .6985E ♦ CO 
0.1694E+00  0.6352F+C0 
-C.2716F+00  3. 4 79  3F ♦ 00 
-0  .701  IF  +00  0.2026E+C0 
-C.1002F+01-0.1654F+C0 
-C.IULF  + 01-0.549CE  + CC 
-0. 1 052 E +01-0. 866 3 F *00 
-C.9062E+00-0. 1074E+CI 
-0.7700F+00-0.1169F+01 


REAL  JP 
0.531 7E+00 
0.5816E+00 
0,8 44  IE + 00 
0. I 145E  + 01 
0 • 1 28  7E  + 01 
0.1 186E+01 
0.9495E+00 
0.7  398E+00 
0.7  135E  + 00 


IMAG  JP 
C.5158E+00 
C.1870E+00 
0.6579E-01 
0.2232E+00 
0.6279E+00 
0. 1077E+01 
C.1383E+01 
0.1495E+01 
0. 1353E+01 
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NhEC  = 
NE  EC  = 


3,  KT * 2,  KR  = 

3*  KT  = 2,  KR  = 


1. 

2, 


SIGHA/(LANBDAl**2x  0.2903E-01 
SIGMA/CL AMBD A) **2=  Q.3845E+00 


VIII.  EXAMPLES  AND  DISCUSSION 
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Some  examples  of  computations  obtained  from  the  computer  program 
are  given  in  the  previous  report  [1],  We  here  give  some  additional  ex- 
amples to  illustrate  the  use  of  the  program. 

As  discussed  in  [1],  the  parameter  a in  the  combined  field  solution 

should  be  taken  0 < a <_  1,  with  a value  of  the  order  of  0.2  normally  giving 

good  results.  To  give  some  additional  information  on  this  choice.  Table 

1 gives  the  RMS  error  in  current  (A)  and  the  backscattering  (a/ira  ) for 

a conducting  sphere  of  radius  a.  The  three  values  ka  = 1.50,  2.75,  and 

4.00  were  chosen  for  ka,  and  the  parameter  a was  varied  from  0.2  to  1.0. 

Also  shown  are  results  for  the  H-field  solution  (a  = 0)  and  for  the  E-field 

2 

solution  (a  = °°)  . The  exact  values  of  cr/ira  are  shown  in  the  last  row,  com- 
puted from  the  spherical  mode  solution  [4].  As  is  evident  from  the  table, 
the  value  a = 0.2  always  gave  good  results,  but  the  error  was  small  for  any 
value  of  a from  0.2  to  1.0. 

2 

Table  1.  RMS  Error  in  current  (A)  and  Backscattering  (cr/ita  ) for  various 
solutions  for  a conducting  sphere  of  radius  a.  H-field  solution  corresponds 
to  a = 0.  E-field  solution  corresponds  to  a » ®. 


RMS  error  in  (A) 

2 

Backscattering  (a/ita  ) 

ka  = 

1.50 

2.75 

4.00 

1.50 

2.75 

4.00 

H-field 

0.0132 

-0.1342 

• 0.0243 

1.076 

0.8702 

0.7274 

a=  0.2 

0.0118 

0.0164 

0.0274 

1.075 

0.8645 

0.7842 

a-  0.4 

0.0117 

0.0187 

0.0362 

1.070 

0.8648 

0.7821 

a=  0.6 

0.0127 

0.0211 

0.0447 

1.065 

0.8647 

0.7912 

o°  0.8 

0.0141 

0.0233 

0.0516 

1.060 

0.8645 

0.8040 

»=  1.0 

0.0154 

0.0254 

0.0571 

1.057 

0.8637 

0.8159 

E-field 

0.0514 

0.5991 

0.1894 

1.043 

0.8246 

0.8586 

Exact 

0 

0 

0 

1.076 

0.8552 

0.7853 

The  problem  of  scattering  by  a conducting  right  circular  cylinder  of 
cross  sectional  radius  a and  height  2a  was  considered  to  illustrate  further 


# 
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the  effect  of  internal  resonance  on  the  various  solutions.  The  exci- 
tation is  a plane  wave  axially  incident  on  the  flat  end  of  (the  cylinder. 

The  first  internal  resonance  varying  as  cos  <J>  occurs  at  ka  = 2.42  approxi- 

2 

mately.  Figure  1 shows  the  normalized  backscattering  cross  section  (cr/X  ) 
in  the  vicinity  of  the  first  internal  resonance.  The  range  of  ka  is  from 
2.37  to  2.47  in  increments  of  0.005.  Note  that  the  internal  resonance  had 
no  observable  effect  on  the  combined  field  solution  with  a - 0.3  (triangles), 
a small  effect  on  the  E-field  solution  (circles) , and  a larger  effect  on  the 
H-field  solution  (squares) . 

To  estimate  the  RMS  error  in  the  current,  the  combined  field  solution 
was  considered  to  be  the  exact  solution.  The  RMS  error  is  then  defined  as 


A = 


|j  - JC|2  ds 


ds 


1/2 


(17) 


where  J is  combined  field  solution  and  the  integration  is  over  the  sur- 
face of  the  scatterer.  Figure  2 shows  A for  the  right  circular  cylinder 
vs.  ka  in  the  vicinity  of  the  first  internal  resonance.  Note  that  the 
greatest  RMS  error  occurs  in  the  E-field  solution  (circles) , with  a some- 
what smaller  error  in  the  H-field  solution  (squares).  However,  since  the 
H-field  eigencurrents  radiate  external  to  the  cylinder,  and  the  E-field 
eigencurrents  do  not  (at  least  theoretically),  the  H-field  backscattering 
cross  section  suffers  a greater  error  at  resonance  than  does  the  E-field 
backscattering  cross  section  (Fig.  1). 

In  the  E-field  solution  a resonance  phenomena  sometimes  occurs  that 
is  unrelated  to  any  internal  resonance.  For  example,  in  the  case  of  a 
sphere  there  was  a small  peak  in  A at  about  ka  = 1.3  (Figure  5,  reference 
[1]).  In  that  case  the  resonance  effect  was  too  small  to  show  up  in  the 
backscattering  cross  section  (Figure  6,  reference  [1]).  This  E-field 
"false  resonance"  effect  was  also  found  to  occur  in  the  case  of  scattering 
by  a flat-back  cone.  The  cone  geometry  is  shown  in  Fig.  3.  The  cone 
geometry  is  shown  in  Fig.  3.  The  excitation  is  a plane  wave  axially 
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2.37  2.42  2.47 


Fig,  1.  Normalized  radar  cross  section  (o/X2)  vs  ka  in  the  vicinity 
of  the  first  cos  $ internal  resonance  (ka  * 2.42)  for  a 
closed  conducting  circular  cylinder  of  radius  a and  height  2a. 
The  H-field  solution  is  shown  by  squares,  the  E-field  solution 
by  circles,  and  the  combined  field  solution  with  a • 0.3  by 
triangles.  The  excitation  is  an  axially  incident  plane  wave. 
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3.27 


Fig.  2.  RMS  error  (A)  in  current  vs.  ka  in  the  vicinity  of  the  first 
cos  4>  internal  resonance  (ka  «*  2. 42)  for  a closed  conducting 
circular  cylinder  of  radius  a and  height  2a,  The  combined 
field  solution  with  a = 0 3 was  assumed  to  be  the  exact  solu- 
tion. The  H-field  solution  is  shown  by  squares  and  the  E- 
field  solution  by  circles.  The  excitation  is  an  axially 
incident  plane  wave. 
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incident  on  the  tip  of  the  cone.  Figure  4 shows  the  normalized  cross 
section  (o/X  ) vs.  frequency  in  the  vicinity  of  an  E-field  false  resonance, 
which  occurred  at  0.9  fQ,  where  fQ  is  the  frequency  at  which  a =*  3X/16. 

Note  that  the  E-field  solution  (circles)  suffers  a considerable  perturba- 
tion at  this  false  resonance,  while  the  H-field  solution  (squares)  and 
the  combined-field  solution  (triangles)  are  unaffected  by  it.  Again  a 
is  taken  to  be  0.3  in  the  combined  field  solution.  It  should  also  be 
noted  that  the  H-field  solution  differs  from  the  combined  field  solution 
by  an  appreciable  amount  in  this  range  of  frequencies. 

Figure  5 shows  the  RMS  error  in  current  for  the  E-field  and  H-field 
solutions  in  the  vicinity  of  the  false  resonance  for  the  flat-back  cone. 
Again  the  combined  field  solution  with  a = 0.3  is  taken  to  be  the  exact 
solution,  with  A computed  according  to  (17).  Note  that  there  is  no 
resonance  effect  apparent  in  the  H-field  solution,  even  though  there  is 
a large  effect  in  the  E-field  solution.  It  is  felt  that  this  false 
resonance  is  a consequence  of  energy  storage  due  to  the  current  being 
approximated  by  triangular  expansion  functions.  If  this  is  a correct 
conjecture,  the  resonance  would  probably  not  be  present  at  the  same 
frequency  if  the  expansion  functions  were  changed,  but  might  occur  at 
a different  frequency. 

In  summary,  we  repeat  that  for  the  combined  field  solution  a should 
be  chosen  in  the  range  0.2  <_  a <_  1.  At  internal  resonances,  the  RMS 
error  in  current  is  greater  for  the  E-field  solution  than  for  the  H-field 
solution.  However,  the  scattering  cross  section  error  is  usually  smaller 
for  the  E-field  solution  than  for  the  H-field  solution.  No  effects  due 
to  internal  resonances  have  ever  been  observed  in  the  combined-field 
solution.  Finally,  the  E-field  solution  sometimes  exhibits  a false 
resonance  effect,  giving  incorrect  solutions  at  some  frequency  which 
does  not  correspond  to  an  internal  resonance. 
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Normalized  radar  cross  section  (a A2)  vs.  frequency  in  the 

vicinity  of  an  E-field  solution  false  resonance  (a  * 3X/16 

at  f ) for  a flat-back  cone  of  radius  a and  90°  cone  angle 
o 

The  H-field  solution  is  shown  by  squares,  the  E-field  solu 
tion  by  circles,  and  the  combined  field  solution  with  a ■ 
by  triangles.  The  excitation  is  a plane  wave  axially  ir.ci 
dent  on  the  cone  tip. 
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Fig.  5.  RMS  error  (A)  in  current  vs.  frequency  in  the  vicinity  of  an 

F,-field  solution  false  resonance  (a  «*  3X/16  at  f ) for  a fl.it- 

o 

back  cone  of  radius  a and  90°  cone  angle.  The  combined  field 
solution  with  a ■ 0.3  waa  assumed  to  be  the  exact  solution. 

The  H-field  solution  is  shown  by  squares  ana  the  F.-field  solu- 
tion by  circles.  The  excitation  is  a plane  wave  axially  inci- 
dent on  the  cone  tip. 
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METRIC  SYSTEM 


•ASE  UNITS: 

Quantity 


length 

mui 

time 

electric  current 
thermodynamic  temperature 
amount  of  aubetance 
luminous  intensity 


SUPPLEMENTARY  UNITS: 


piane  angle 
•olid  angle 

DERIVED  UNITS 

Acceleration 

ect'vitv  (of  a radioactive  source) 
angular  acceleration 
angular  velocity 
tree 

dnrSI'v 

electric  capacitance 

electrical  conductance 

elertnc  field  strength 

electric  inductance 

electric  potential  difference 

electric  resistance 

electromotive  force 

energy 

entropy 

force 

frequency 

illuminance 

luminance 

luminous  flux 

magnetic  field  strength 

magnetic  flux 

magnetic  flux  density 

magnetomotive  force 

power 

pressure 

i|Ui»’!tit>  >1  electricity 
quantity  of  heat 
radiant  intenaitv 
sp**i  dir  heat 
streaa 

thermal  uinducfivity 
velocity 

viscosity  dynamic 

viscosity,  kinematic 

voltage 

volume 

wavenumber 

work 


Unit 


metre 

kilogram 

second 

ampere 

kelvm 

mole 

candela 


radian 

steradian 


metre  per  second  squared 

disintegration  per  second 

radian  per  second  squared 

radian  per  second 

square  metre 

kilogram  per  cubic  metre 

farad 

siemens 

volt  per  metre 

henry 

volt 

ohm 

volt 

■oule 

toule  per  kelvin 

newton 

hertz 

lux 

candela  per  square  metre 
lumen 

ampere  per  metre 

weber 

tesla 

ampere 

watt 

postal 

coulomb 

loule 

watt  per  steradian 

mule  ; nr  ktlogram-kelvin 

pascal 

watt  per  metre  kelvin 
metre  per  second 
pascal-second 
square  metre  per  second 
volt 

cubic  metre 
reciprocal  metre 
touie 


SI  Symbol 
m 

k« 

s 

A 

K 

mol 

cd 


rad 

tr 


F 

S 

H 

V 

V 

I 

N 

Hz 

lx 

Iro 

Wb 

T 

A 

W 

Pa 

C 

I 


Pa 


V 


i 


Formula 


ml* 

(disintegration)/* 

rad/s 

rad/s 

m 

kg/m 

A -*N 

AN 

Vim 

V-a/A 

W/A 

V/A 

W/A 

N-m 

VK 

kg- m/a 

(cyclef/e 

bn/m 

cd/m 

akr 

Aim 

V-s 

Wbrtn 

ji» 

N/m 

A-a 

N-m 

W/ar 

J/kg-K 

N/m 

W/m-K 

ml* 

Pet 

ml* 

W/A 

m 

(waveVm 

N-m 


SI  PREFIXES: 


Multiple  noun  Factors 

Prefix 

SI  Sym 

1 000  000  IKK)  000 

10" 

ten 

T 

<; 

1 000  IKK)  000  ~ 

10’ 

til* 

1 IH  0 IKK)  • 

10* 

tnaga 

M 

t 000  =- 

10’ 

kite 

k 

too  - 

10* 

hacto* 

h 

10  - 

10' 

deka’ 

da 

0 t 

10  1 

dad* 

d 

0 01  - 

10  1 

naotl* 

c 

.1001 

HI"' 

mllll 

m 

ii  iNKi  mu 

10* 

micro 

M 

0 fHHI'HKI  (K‘1 

itl  * 

nano 

n 

0 000  IMHl  IKK)  (KU 
0 000  OOtl  IKK)  (MKl  mu 

10  " 
10  " 

ptco 
mm  to 

r 

IKK1  (MM)  000  IKK!  IKK)  fKU 

10  " 

aflo 

• 

‘ T j be  avoided  wnere  possible 
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